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Abstract. We investigate the flow and jamming through small apertures of a column 
of granular disks via a pseudo-dynamic model. We focus on the effect that the 
preparation of the granular assembly has on the size of the avalanches obtained. 
Ensembles of packings with different mean packing fractions are created by tapping the 
system at different intensities. Surprisingly, packing fraction is not a good indicator 
of the ability of the deposit to jam a given orifice. Different mean avalanche sizes 
are obtained for deposits with the same mean packing fraction that were prepared 
with very different tap intensities. It has been speculated that the number and size 
of arches in the bulk of the granular column should be correlated with the ability of 
the system to jam a small opening. We show that this correlation, if exists, is rather 
poor. A comparison between bulk arches and jamming arches (i.e., arches that block 
the opening) reveals that the aperture imposes a lower cut-off on the horizontal span 
of the arches which is greater than the actual size of the opening. This is related to 
the fact that blocking arches have to have the appropriate orientation to fit the gap 
between two piles of grains resting at each side of the aperture. 
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1. Introduction 

The flow of granular materials through apertures is commonplace in a variety of 
industrial applications. Studies in this respect can be separated into two major areas: 
(1) continuous flow, and (2) jamming. Continuous granular flow is observed for dry non- 
cohesive materials if the size of the aperture is large (typically above five grain diameters 
for spherical particles). Jamming is observed whenever the opening is smaller, which 
requires the input of external perturbations in order to restart the flow. The cause of 
jamming is the formation of a blocking arch. We distinguish two type of arches: blocking 
arches (or jamming arches) and bulk arches. Blocking arches form at the orifice during 
discharge and prevent the flow. Bulk arches form at any place inside the packing during 
the dynamical process that leads the grains to reach mechanical equilibrium. Arches 
are set of particles that are mutually stable. The removal of any particle in the arch 
leads to the destabilization of the others. 

A number of studies have considered the jamming of an aperture during the 
discharge of grains. These include experimental studies on two-dimensional (2D) 
hoppers using circular grains [Tl[2] and three-dimensional (3D) silos using spherical and 
non-spherical particles [3111], numerical simulations using discrete element methods in 
2D [ZllH], experiments with quasi-2D silos and spherical grains |5], 3D vibrated silos [6] 
and experiments with tilted [9] and wedge-shaped hoppers [lO]. Also, the properties of 
blocking arches and bulk arches have been considered in the past [Tl fTTHTS] . However, 
none of these studies have considered the effect of the preparation of the granular column 
prior to the discharge. It is known that the number and size of arches inside a granular 
assembly are dependent on the packing fraction. Therefore, it is expected that the 
jamming of columns prepared at different packing fractions may occur with different 
probability. A related issue is the question of to what extent the arches formed in the 
bulk of the system are comparable with the arches that effectively block the aperture 
during drainage. 

In this paper we use a 2D pseudo-dynamic simulation scheme previously developed 
by Manna and Khakhar [16l[T7] to study the effect of the initial packing fraction on the 
jamming probability and the correlation between bulk arches and blocking arches. The 
jamming probability is directly connected with the mean size of the avalanches [5]. An 
avalanche is the flow of grains that occurs between the initiation of the discharge and 
the arrest of the flow due to the formation of a blocking arch [3]. We will show that 
there is a strong dependence of the size of the avalanches with the packing fraction. 
However, there is not a monotonic relation between these two quantities. We also flnd 
that there is a poor correlation between the size of arches in the bulk and the size of 
the avalanches. A comparison between bulk arches and jamming arches reveals that 
the aperture not only imposes a cut-off on arches of horizontal span below the opening 
size, but also prevents the formation of some blocking arches that, in principle, are wide 
enough to induce jamming. 
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2. The pseudo-dynamic algorithm 

Our simulations are based on an algorithm for inelastic massless hard disks designed by 
Manna and Khakhar [T6l[T7] . This is a pseudo-dynamics that consists in small falls and 
rolls of the grains until they come to rest by contacting other particles or the system 
boundaries. We use a container formed by a fiat base and two fiat vertical walls. No 
periodic boundary conditions are applied. 

The deposition algorithm consists in choosing a disk in the system and allowing a 
free fall of length 6 if the disk has no supporting contacts, or a roll of arc-length 5 over 
its supporting disk if the disk has one single supporting contact [T6] - [T8] . Disks with two 
supporting contacts are considered stable and left in their positions. If in the course of 
a fall of length 6 a disk collides with another disk (or the base), the falling disk is put 
just in contact and this contact is defined as its first supporting contact. Analogously, 
if in the course of a roll of length 6 a disk collides with another disk (or a wall), the 
rolling disk is put just in contact. If the first supporting contact and the second contact 
are such that the disk is in a stable position, the second contact is defined as the second 
supporting contact; otherwise, the lowest of the two contacting particle is taken as the 
first supporting contact of the rolling disk and the second supporting contact is left 
undefined. If, during a roll, a particle reaches a lower position than the supporting 
particle over which it is rolling, its first supporting contact is left undefined (in this way 
the particle will fall vertically in the next step instead of rolling underneath the first 
contact). A moving disk can change the stability state of other disks supported by it, 
therefore, this information is updated after each move. The deposition is over once each 
particle in the system has both supporting contacts defined or is in contact with the 
base (particles at the base are supported by a single contact). Then, the coordinates 
of the centers of the disks and the corresponding labels of the two supporting particles, 
wall, or base, are saved for analysis. 

An important point in these simulations is the effect that the parameter 6 has 
in the results since particles do not move simultaneously but one at a time. One 
might expect that in the limit 6 ^ we should recover a fairly "realistic" dynamics 
for fully inelastic non-slipping disk dragged downwards at constant velocity. This 
should represent particles deposited in a viscous medium or carried by a conveyor belt. 
Although this dynamics contrasts with the dynamics of dry granulates, experiments on 
the jamming of fully submerged grains in gels [I9] have shown remarkable similarities 
with the more widely available data on dry systems. We chose 6 = 0.0062d (with d 
the particle diameter) since we have observed that for smaller values of 6 results are 
indistinguishable from those obtained here |18j . 

The pseudo-dynamics approach has been chosen in view of the low CPU time 
demanded by this scheme. In this study, we need to prepare a large number of samples 
(through tapping) with a given packing fraction and then trigger discharges for different 
aperture sizes. More realistic simulations such as granular dynamics (or discrete element 
method) would require a much higher computational effort. In spite of the simplifications 
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Figure 1. (Color online). Sample eonfigurations of the granular column for the 
steady state corresponding to F — 0.39 before (a) and after (b) discharging through 
and opening of width D = 2.75. The red segments indicate the arches in the system. 



of the pseudo-dynamics, it has been shown that results on tapping agree quahtatively 
with granular dynamics simulations |21] . 

3. Initial packings 

In order to study the effect of packing fraction on the jamming of the flow through an 
aperture, we first need to prepare packings at reproducible packing fractions. To achieve 
this, several techniques can be applied. For example, the sequential deposition of grains 
submerged in a viscous liquid yield reproducible packing fractions that can be tuned by 
changing the friction coefficient of the particles or the density mismatch |22]. We have 
chosen another well known technique to generate reproducible ensembles of packings: 
tapping. Nowak et al. have shown that an appropriate tapping protocol can lead to 
reproducible states in the sense that an ensemble of configurations with well defined 
mean packing fraction is recovered if the same protocol is followed irrespective of the 
initial state [23]. This has been more carefully discussed by Ribiere et al. [20]. Dijksman 
et al. showed how different states can be obtained not only by changing the tap intensity 
but also by changing the tap duration [23]. A similar effect was investigated by Pica 
Ciamarra et al. in submerged samples where a fluid pulse is used to excite the granular 
column [25]. Hence, we use a simulated tapping protocol (see below) to generate sets of 
initial configurations that have well defined mean packing fractions. 

The simulations are carried out in a rectangular box of width 24.78rf containing 1500 
equal-sized disks of diameter d. Initially, disks are placed at random in the simulation 



Effect of packing fra 



0.84 
0.82 
0.8 
'0.78 
0.76 
0.74 
0.72, 







oooo 



0.5 



1.5 



Figure 2. Mean packing fraction in the steady state of the tapping protocol as a 
function of the tap intensity F. 



box (with no overlaps) and deposited using the pseudo-dynamic algorithm. Once all 
the grains come to rest, the system is expanded in the vertical direction and randomly 
shaken to simulate a vertical tap. Then, a new deposition cycle begins. After many 
taps of given amplitude, the system achieves a steady state where all characterizing 
parameters fluctuate around equilibrium values independently of the previous history 
of the granular bed. The existence of such "equilibrium" states has been previously 
reported in experiments |20j . 

The tapping of the system is simulated by multiplying the vertical coordinate of each 
particle by a factor A (with ^4 > 1). Then, the particles are subjected to several (about 
20) Monte Carlo loops where positions are changed by displacing particles a random 
length Ar uniformly distributed in the range < Ar < A — 1. New configurations 
that correspond to overlaps are rejected. This disordering phase is crucial to avoid 
particles falling back again into the same positions. Moreover, the upper limit for Ar 
(i.e. ^ — 1) is deliberately chosen so that a larger tap promotes larger random changes in 
the particle positions. The expansion amplitude A ranges from 1.03 up to 3.0. Following 
Refs. |211[26] we quantify the tap intensity by the parameter F = y/A — 1. For each 
value of F studied, 10^ taps are carried out for equilibration followed by 5 x 10'^ taps for 
production. 500 deposited configurations are stored which are obtained by saving every 
10 taps during the production run after equilibration. These deposits will be used later 
as initial conditions for the discharge and flow through an opening. 

The deposited configurations are analyzed in search of bulk arches. We first identify 
all mutually stable particles — which we define as directly connected — and then we find 
the arches as chains of connected particles. Two disks A and B are mutually stable if A is 
the left supporting particle of B and B is the right supporting particle of A, or viceversa. 
We measure the total number of arches, arch size distribution n{k), and the horizontal 
span distribution of the arches nk{x). The latter is the probability density of finding an 
arch consisting of k disks with horizontal span between x and x + dx. The horizontal 
span (or lateral extension) is defined as the projection onto the horizontal axis of the 
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Figure 3. (Color online). Number of arches per particle (red circles) and mean size 
of the arches in terms of the number of grains involved in an arch (black squares) as a 
function of F. 



segment that joins the centers of the right-end disk and the left-end disk in the arch. 
In Fig. [Dj^a), we show an example of a deposited configuration with arches indicated 
by segments (for a description of Fig. [U^b) see next section). Notice that the pseudo- 
dynamics mimics the behavior of disks that roll without slipping. This corresponds to a 
system with infinite static friction which is expected to yield a large number and variety 
of arches. The arch structure of frictionless systems may differ significantly from the one 
seen in Fig. [H^a). However, simulations with realistic discrete element methods with 
finite friction yield similar structures |13] . 

In Fig. [21 we present the steady state mean packing fraction, 0, of our granular 
deposits as a function of F. There exists a rather sharp decrease of as the tapping 
intensity is increased followed by a minimum and a very smooth increase. The sharp 
drop of is associated to a discontinuous order- disorder transition previously reported 
for this model [18] and also observed in granular dynamics of polygonal grains [27]. 
The appearance of the minimum packing fraction as a function of tap intensity has 
been reported for several models (including a frustrated lattice gas model [28], a Monte 
Carlo type deposition [21] and a realistic discrete element method simulation [21]) and 
in experiments of tapping with a quasi-2D system [29l[30]. For the model we use in 
this paper, the minimum has been shown to exist even if bidisperse systems are 
considered [3T] . 

The minimum in is caused to a large extent by the formation of arches [21] . Figure 
|3] shows the number and mean size of arches found in the system as a function of F. 
When F is increased considerably (above 0.7), every tap expands the assembly in such a 
way that particles get well apart from each other. During deposition, particles will reach 
the free surface of the bed almost sequentially (one at a time), reducing the chances of 
mutual stabilization. Therefore, arches are less probable to form as F increases and so 
must grow since fewer voids get trapped. Indeed, we see in Fig. [3] that the number 
and size of arches decrease at large F (> 0.7) for increasing tap intensities. Eventually, 
for very large F, no arches are formed after each tap and will reach a limiting value. 
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At lower F (0.4 ^ T < 0.7), the free volume injection due to a tap creates very narrow 
gaps between particles. For a given arch to grow by the insertion of a new particle, it 
is necessary to create a gap between two particles in the existing arch where the new 
particle can fit in. This explains why increasing F will promote the formation of larger 
arches and reduce (j) in this regime. For very low tapping intensities (F < 0.3), we find 
a rather constant (p. However, the number and size of the arches decrease with F (see 
Fig. [3]). This would imply that a maximum in the packing fraction should be observed at 
such light tapping. This maximum has been recently reported in other models [271132]. 
but is not present in our pseudo-dynamics. Notice that the number and size of the 
arches give only a rough indication for the free volume in the sample since the actual 
shape of the arches will also be important. 

4. Flow and jamming 

For each deposit generated as described in the previous section, we trigger a discharge by 
opening an aperture of width D relative to the diameter d of the disk in the center of the 
base of the containing box. Grains will flow out of the box following the pseudo-dynamics 
until a blocking arch forms or until the entire system is discharged (with the exception of 
two piles resting on each side of the aperture). During the dynamics, disks that reach the 
bottom and whose centers lie on the interval that defines the opening will fall vertically 
(even if the surface of the disk touches the edge of the aperture). This prevents the 
formation of arches with end disks sustained by the vertical edge of the orifice. Although 
such arrangements happen in real experiments, they are uncommon [Tl]. After each 
discharge, we record the size of the avalanche (i.e. the number of grains flowed out) 
and the final arrangement of the grains left in the box. Averages are taken over 500 
discharges for each value of F used to prepare the initial packings. 

One single discharge attempt is carried out for each initial deposit. This allows us 
to assure that the initial preparation of the pack belongs to the ensemble of deposits 
corresponding to the steady state of the particular tap intensity chosen. In many 
experiments and most industrial applications, discharges are triggered one after another 
from the same deposit without preparing the system in the initial condition again [3H6] . 
However, some experiments do fill the container anew before each discharge [ni2|l33]. 
We can see in Fig. [11(b) an image of the system after a discharge that resulted in a jam. 
It is clear that the structure of the packing is greatly affected by the partial discharge 
in our simulations. Therefore, these final structures are not used for new discharges. 

In Fig. [H we plot the avalanche size distribution p{s) for a few values of F and 
D = 2.25. We obtain this by counting the number of grains s that flow in each of the 
500 discharges corresponding to each initial deposit generated for each F. An exponential 
tail in p{s) has been observed in several previous studies, both two-dimensional [2l[5l[71[8] 
and three-dimensional [3l[l]. Manna and Herrmann [3l], using the same model, found 
avalanches with a power law distribution. Notice however, that in Ref. [M] the authors 
trigger one avalanche after the other by simply removing a grain of the blocking 
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Figure 4. (Color online). Avalanche size distribution p{s) for several values of F for 
D = 2.25. (a) Semilog plot, (b) log-log plot. 



arch. This minute perturbation to trigger avalanches may induce strong correlations 
between successive discharges in contrast with the strong rearrangements induced in 
most experiments. Based in our hmited number of discharges, we are unable to assert if 
an exponential or a power-law decay is at play in our simulations (see log-lin and log-log 
plots in Fig. H]). Although we report p{s) up to s = 100, larger avalanches of up to 
1000 disks are observed, whose statistics can be largely affected by the finite size of the 
system {N = 1500). 

The mean avalanche size, (s), as a function of F is shown in Fig. |5] for various 
aperture sizes. As it can be expected, (s) increases if D increases. As we can see, for 
small apertures, (s) grows monotonically as F increases. However, for D > 2.0, the 
mean avalanche size presents a local maximum and a local minimum as a function of F. 

It has been speculated [TT1|35] that the size of the avalanches can be connected 
with the arches inside the granular deposit. Although arches in the initial configuration 
are not dragged to the aperture during flow since arches actually break and form all 
the time in the process, it is believed that the ability of the system to form arches in 
the initial deposit is connected with the ability to form blocking arches during flow. 
Indeed, the features observed in Fig. |5]are somewhat correlated with the number and 
size of the bulk arches. As F is increased from the lower values, the size of the arches 
remains initially rather constant (whereas the number of arches decreases. Fig. |3]). This 
reduction in the number of arches leads to a smaller jamming probability and a rapid 
increase in the size of the avalanches (see Fig. |5]). This regime ends when the sharp drop 
in ends (F ^ 0.4). For larger F and up to the packings with minimum packing fractions 
(i.e. 0.4 < F < 0.7), the number and size of the arches increase. As a consequence, (s) 
decreases due to the increased likelihood of jamming. Finally, for F > 0.7, the number 
and size of the arches fall and a corresponding increase of the avalanche sizes is observed. 
From these observations we can assert that the number and size of arches in a given 
packing give an overall indication of the chances that the system will jam if it is left to 
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Figure 5. (Color online). Mean avalanche size (s) as a function of F for several 
sizes of the aperture D. Notice that (s) is affected to a large extent by the rare large 
avalanches not reported in Fig. 31 



flow through a small aperture. 

Notwithstanding the previous analysis, the size of the avalanches is not only 
dependent on the preparation protocol — defined in our case by F — but also on the 
actual size of the outlet imposed. For example, as we mentioned, for D = 2.0 the mean 
avalanche size does not present the maximum and minimum suggested by the number 
and size of arches. In order to take into account the effect of the size of the aperture 
we measure the probability Parch of finding an arch in the bulk of the deposits wide 
enough to block a given aperture D. We measure the horizontal span of each arch as 
the projection on the horizontal axis of the segment that joins the centers of the end 
particles of the arch. An arch of span x can jam an opening of width x + d (with d the 
diameter of a grain). In this analysis we include the grains that do not form arches, 
which can jam any orifice with D < 1. In Fig. [6] we plot Parch(-D) as a function of F 
and compare with the corresponding (s). Overall, the probability of finding an arch 
wide enough to block an aperture of size D decreases with D in correspondence with 
the overall increase of (s). However, for a given D, the dependence with F does not 
show a clear anti-correlation between the probability and the mean avalanche size. This 
implies that the bulk arches can give only a rough indication of the eventual size of the 
avalanches that would discharge if an aperture is opened. 

It is important to mention that we have always opened the aperture at the center 
of the bottom of the container. The ordering observed for F < 0.5 (see Fig. [1]) suggests 
that an effect related to the relative position of the aperture and the first layer of grains 
may be expected. The main effect would be related to the fact that particles just at 
the edges of the orifice do not flow in the pseudo-dynamics and these will produce a 
reduced effective aperture. If all discharges start from an initial packing so ordered that 
grains at the first layer always sit on the same horizontal positions, the mean avalanche 
size would depend on the horizontal position of the aperture. When this happens, 
one observes oscillations in the avalanche size as a function of the aperture size |33j . 






Figure 6. (Color online). Mean avalanche size (s) as a function of T and probability 
Paich{D) of finding an arch of horizontal span _D — 1 for several apertures: (a) D = 2.0, 
(b) D = 2.25, (c) D = 2.5, (d) D = 2.75. 



However, this does not happen in our simulations as we can see in Fig. |5l This effect is 
observed only for highly ordered structures obtained with frictionless particles [33] . Our 
particles model non-slipping grains and the small deviations in position of the disks of 
the first layer with respect to a truly crystalline structure are sufficiently large to mask 
any systematic effect due to ordering that may require a detailed study on the position 
chosen for the aperture. 

5. Effect of packing fraction 

It is generally believed that packing fraction is a good parameter to characterize many 
properties of a granular bed [36l[37]. Results obtained for deposits prepared at a given 
are not necessarily general and must be repeated for different packing fractions. 
However, this does not mean that packing fraction is the only or the main factor that 
can affect the results. In Table [1] we present part of the data of Fig. [5] but ordered by 
mean packing fractions (p corresponding to the steady state obtained for each given F 
(see Fig. |2]). As we can see, the mean avalanche size depends on (;/> in a non trivial way. 
At the highest 0, obtained by light tapping, the mean avalanche size can range from a 
few grains to a hundred grains depending on the value of F used to create the packings. 
The larger the aperture D, the wider the range in (s) within this regime where the 
system is rather ordered but the number and size of the arches fall with increasing F 
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Table 1. Mean avalanche size (s) for different tap intensities F that yield similar 
packing fractions 0. 

(see Fig. [3]). For low deposits with similar packing fractions but created with low 
and high tap intensities display different values of (s) for any given D. 

As we can see, deposits with the same </> can present different values of (s). 
Therefore, the steady state ensembles of packings with equal obtained by tapping 
may behave differently. This has been pointed out in Refs. [29l[30] where steady states 
with the same but bearing different stresses were obtained through tapping. In our 
simulations, forces are not calculated and therefore the stress tensor cannot be obtained. 
However, a clear difference in the response of the granular columns with same is 
observed in the sense that avalanches are, in average, of different size. 

6. Connection between bulk arches and jamming arches 

Although we have shown in the previous section that arches found in the bulk of 
the granular packing give a rough indication as to whether the system would be 
more or less likely to jam, arches actually formed at the aperture during discharge 
are different. Detailed studies of such blocking arches have been reported for two- 
dimensional experimental setups [131I38] . Since we have access to both bulk and blocking 
arches in our simulations, we compare a few properties and discuss on the implications 
for the jamming probability. 

In Fig. [7] we compare the arch size distribution nifi) for the bulk arches and for the 
jamming arches for different values of D. n{k) is the probability of finding an arch of k 
grains (A; > 2). For bulk arches, n{k) is calculated as the number of arches of k disks 
found in all initial packings divided by the total number of arches (i.e., summing for all 
k > 2). For jamming arches, n{k) is calculated as the number of discharges that led to 
a blocking arch of k disks divided by the total number of discharges (discharges that 
ended without producing a jam were not considered). The plot is presented as a function 
of k — D since this produces a collapse of the curves by subtracting a quantity (D) 
proportional to the lower cut-off imposed by the orifice on the arch sizes (see Ref. |14j). 
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Figure 7. (Color online). Distribution n{k) of arch sizes for F = 0.387. The black 
symbols correspond to arches found in the bulk prior to the discharge, the color data 
correspond to the blocking arches for different values of D as indicated in the legend. 
Notice that the horizontal axis is shifted by the size of the aperture D for each set of 
data. 



The fact that all the normalized histograms fall on the same curve indicates that the 
nature of the arches that jam the aperture is the same for small and big orifices. In 
general, an exponential tail is observed and a cut off for small arches is imposed by 
the orifice. The large arches generally form upstream supported by two stationary piles 
resting on both sides of the opening. It is clear from Fig. [7] that blocking arches tend 
to be larger than bulk arches (even after the correction due to the cut-off imposed by 
D). We believe this is due to the fact that many small arches that are stable in the 
bulk thanks to the many neighbors cannot accommodate in the conical shaped funnel 
created by the two stationary piles. This is best demonstrated in the next paragraph. 

The horizontal span of the arches for a given number of grains k is shown in Fig. |H1 
We include data for the bulk arches found in the initial deposits and for the jamming 
arches found for an aperture of width D = 2.5. As we can see, for a given k, arches are 
more likely to be wider in the case of jamming arches as compared to bulk arches. This 
is to be expected for small k since small arches with small span might not be able to 
jam the aperture. However, even arches with bigger k are biased in the distribution of 
blocking arches. This is due to the two piles formed at each side of the orifice. Blocking 
arches that are wider than the aperture D must span this funnel. Arches of k disks with 
horizontal span sufficient to jam the orifice might still be unable to span the funnel (see 
the inset in Fig. E^c)). This results in blocking arches generally wider, for a given k, 
than the corresponding bulk arches. Notice that arches of horizontal span x < D — 1 
may also jam the orifice (see Fig. [Ht^a)) due to the grains sitting at the edges of the 
aperture that reduce the effective size of the opening. This effect has been studied in 
more detail by Pournin et al. |33j . 
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Figure 8. (Color online). Horizontal span distribution nk{x) for arches formed by 
different number of grains fc at F = 0.707. (a) fc = 2, (b) fc = 3, (c) fc = 4, (d) fc = 5. 
The red dashed lines correspond to blocking arches for an aperture D = 2.5, whereas 
the black solid lines correspond to arches found in the bulk of the initial deposits. The 
inset of panel (c) shows two schematic arches of fc = 4 which are wide enough to jam 
the orifice but only one of them fits in the gap left by the two piles at rest on each side 
of the aperture. 



7. Conclusions 

We have considered granular avalanches discharged though small apertures at the 
bottom of a container by using a 2D pseudo-dynamic model. We have focused on the 
effect of the packing fraction of the granular deposit prior to the avalanche discharge. 
The results indicate that the initial packing fraction has an important effect on the mean 
avalanche size for a given opening size. However, similar (s) can be obtained for packings 
with very different (j). Most importantly, very different values of (s) can correspond to 
initial packings with the same packing fractions that were prepared by using different 
tap intensities. It is important to note that these results are obtained not for single 
packings but for ensembles of deposits representative of steady states corresponding to 
a particular tap intensity. 

Our main conclusion is that packing fraction is not a good macroscopic parameter 
to predict the size of the avalanches that would flow through a given aperture. It seems 
that further information is necessary. Although this information is expected to reside 
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in the size and number of arches, we have seen that the correlation of these with (s) is 
not consistent for all openings D. It seems that is not possible to predict the jamming 
probability of a granular column as it flows through a small aperture based on a few 
global properties of the initial deposits. 

A side result from our study is that blocking arches are generally wider than the 
arches found in bulk. This is not only due to the fact that the aperture imposes a lower 
cut-off for the possible jamming arches, but also to the fact that even arches formed in 
the bulk which are wide enough to block the outlet have shapes not compatible with the 
conical boundary effectively created by the two piles of stationary grains at the sides of 
the aperture. 
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